Study of Multiquarks Systems in a Chiral Quark Model 
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We discuss the stability of multiquark systems within the recent model of Glozman 
et a], where the chromomagnetic hyperfine interaction is replaced by pseudoscalar- 
meson exchange contributions. In this model the {u, d) diquark 5 = / = is 
strongly bound (more than for the chromomagnetic interaction) and this leads to 
a bound heavy tetraquark QQqq where the light quarks are in the S = / = 
configuration. We study the stability of other multiquark systems as well. 



1 Introduction 

The interest for the possible existence of multiquark hadrons has been raised 
twenty years ago by Jaffe, wJio suggested that states of two quarks - two 
antiquarksB and of six quarks □ could be bound. 

In the following years this problem has been studied within a large vari- 
ety of models. Some earlier studies in MIT bag indicated the presence of a 
dense spectriim of tetraquark states in the light sectorld (and more generally of 
multiquarkslj) . Later on, tetraquark systems have been examined in potential 
models oO'Bfl and flux tube modelsB. Weinstein and Isgur showedQ that there 
are only a few weakly bound states of resonant meson-meson structure in the 
light {u, d, s) sector. Usually the ground state of these systems lies closely 
above the lowest (qq) + {qq) threshold. The other extreme result, asijompared 
to MIT bag model, was obtained by Carlson and Pandharipande H in their 
flux-tube model with quarks of equal masses, where no bound state was found. 

Altogether the theoretical predictions about the existence of these bound 
states are still unclear. 

As far as the experimental situation is concerned, there are some candi- 
dates for non-qq states, but data are not yet conclusive. Furthermore, even 
if a resonance would be clearly identified as an exotic (not qq or qqq) state, 
nevertheless a careful phenomenological study would be necessary in order to 
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understand its internal structure among many different possibilities (multi- 
quarks, hybrids, glueballs,...). 

The description of these candidates and of the phenomenological properties 
which permit a distinction between diffemal exotics is beyond the purposes of 
this proceeding arui-we refer to Ref.s EjO and the last issue of Review of 
Particle PropertiesEJ for further details. 

Anyway, a well-established theoretical result is that systems, jwith a larger 
mass difference among their components are more easily boundu'cl. For exam- 
ple, for a system of two heavy quarks and two light antiquarks QQqq {Q = c 
or 6, g = M, d or s) stability can be achieved-without spin-spin interaction, 
provided the mass ratio m{Q)/m{q) is larger El than about 15, which means 
that Q must be a 6-quark. 

Recently, the interest for multiquark systems containing charm quarks has 
grown, considering that new experiments are being planned at Fermilab and 
CERN|_tQ_smrch for new hadrons and in particular for doubly charmed tetra- 

quarksyye. 

In this context, we have carried outll3 a study of the. QQnq system in the 
framework of the chiral quark model of Glozman et al. tSilj. This model is 
somehow quite "extreme" because it includes meson-exchange forces between 
quarks and entirely neglects the chromomagnetic interaction. However it per- 
mits a very good description of the baryon spectrum, it is thus worth testing 
it in further predictions. 

Considering that in this model light-quark mesons are "quasiparticles" 
with a spectrum which must be assumed and cannot be evaluated directly (this 
is of course rather an unpleasant feature of the model), multiquark systems 
represent an obvious testing ground for it. Other possibleptests could come 
from a careful analysis of the spectrum of charmed baryonstZl. 

2 The Glozman model 

Before presenting our results about multiquarks, let us briefly consider the 
Glozman et al. model and compare it with other potential models used in 
hadron spectroscopy. 

In a general Hamiltonian, which would approximate the low energy limit 
of QCD, one can introduce both a chromomagnetic interaction and a meson- 
exchange contribution, obtaining an explicit form as 



i<j i<j 

where rrii is the constituent mass of the quark located at r^; rij = |r*j — r^l 
denotes the interquark distance; al, AJ, Af are the spin, colour and flavour op- 
erators, respectively. Spin-orbit and tensor components may supplement the 
above spin-spin forces for studying orbital excitations (they give no contribu- 
tion for L = systems). The potential in (1) has three parts containing the 
confining, the chromomagnetic and the meson-exchange contribution. 

Usually, the confining term Vi^onf is assumed to include a Coulomb plus a 
linear term, 

Vconl = --+ br + c. (2) 
r 

In theJpUowing, we shall either use the very weak linear potential of Glozman 
et al. E3 corresponding to 

(Ci) a==c = 0, and 5 = 0.01839 GeV^ (3) 

or the more conventional choice 

(C2) a 0.5203, 6 = 0.1857 GeV^ c = -0.9135 GeV, (4) 

which has already been applied to the study of tetraquarks by Silvestre-Brac 
and SemaylS. 

The term J2i<j K'^'j ^i'^j ^(^ij) is the chromomagnetic analogue of the 
Breit-Fermi term of QED. In the interaction between a quark and an antiquark 
(as e.g. in a meson) one finds A5-A2 = —16/3. Then a positive Vg shifts each 
vector meson above its pseudoscalar partner, for instance D* > D in the charm 
sector. For baryons, where A^-A^ = —8/3 for each quark pair, such a positive 
Vg pushes the spin 3/2 ground states up, and the spin 1/2 down, for instance 
A > A'^. For the radial shape, as an example, we mention 

V,= "/-P-/^ (5) 

which was used in Ref. i, with the same value of a as in Eq. (4) and d = 
0.454 GeV-i. 

Finally, the last term of H corresponds to meson exchange, and an explicit 
sum over F is understood. If the system contains light quarks only (as in 
Ref.sEaO), the sum over F runs from to 8, i.e. over the members of the 
J^'-^ = ^ nonet, which represent the Goldstone bosons of the spontaneusly 
broken SU{'i)A symmetry (1 — 3 — > tt, 4—7 ^ K, 8 ^ rj and — > 77'). If a heavy 



flavour is incorporated, a phcnomenological extension from SU(3)f to SU(4)f 
would further extend the sum to F = 9 — 12 corresponding to a iJ-exchange, 
F = 13 — 14 to a D^-exchange and F = 15 to an 77c-exchange (of course in this 
case the interpretation as Goldstone bosons is not really possible). Similar 
terms should then be introduced if one includes the beauty sector as well. The 
radial form of 14^ 7^ is derived from the usual pion-exchange potential which 
contains a long-range part and a short-range one 



V- . ^ . ^ .9' 1 



■^,exp(-^_^^^(3)^^^^^- 



(6) 



where /i is the pion mass. 

When constructing AW forces from meson exchanges, one usually disre- 
gards the short-range term in Eq. (6), for it is hidden by the hard core, and 
anyhow the poteiitial in that region is parameterized empirically. For example. 



when Tornqvistkil, Manohar and Wisela or Ericson and KarlB considered pion 
exchange in multiquark states, they used the Yukawa term exp(— /ir)/r acting 
between two well-separated quark clusters. Weber et ai. C3, in their model 
with hyperfine plus pion-exchange interaction, studied both the cases with or 
without delta-term, showing that good results could be obtained without it in 
Jbati^on spectroscopy. Thus the relevance of an ad-hoc regularized delta-term 
E3o in the study of baryon spectroscopy is somehow surprising. Nevertheless, 
this ansatz permits to give a good description of the baryon spectrum and, 
in particular, allows one to solve the problem of the ordering of the lowest 
parity-odd and parity-even states of N , A and S rcsoiiances, which did not 
find a solution in conventional chromomagnetic modelslHj. 

Incidentally, this result is not completely unexpected, in fact Buchman et 
al. had shown, already some years ago, that this term is essential for obtaining 
a good description of magnetic momentsEj. 

For the sake ofxompleteness we report here the explicit form of the regu- 
larized delta-terniEj 

V, = Q{r ro)/.2^^P^^ - ^ exp(-.2(r - ro)'), (7) 

where tq = 2.18 GeV-\ e = 0.573 GeV, and ^ = 0.139 GeV for tt, 0.547 GeV 
for 1] and 0.958 GeV for 77'. Furthermore, the Yukawa- type part is cut off for 
r < vq. This smearing should account for the fact that both the pseudoscalar 
mesons and the constituent quarks have a finite size and that boson fields 
cannot be described by a linear equation near their source. |— , 
The explicit form of the Hamiltaman which extends the results of Ref. Il3 
from SU(3) to SU(4) is given in Ref. ta. It includes also the exchange of D, 



and mesons. However, considering that Httle uu or dd mixing is expected 
in r]c, this contribution can be neglected when considering systems with none 
or one charm quark. Moreover when the meson mass fi reaches values of a few 
GeV as for D or r/c the two terms in_Eq. (6) basically cancel each otheii.and one 
recovers practically the SU(3) formE£l. This is in agreement with Ref.Ej where 
it has been explicitly shown that the dominant contribution to the Sc and E* 
masses is due to meson exchange between light quarks and the contribution of 
the matrix elements with the D (Ds) and D* (D*) quantum numbers (which 
are evaluated phenomenologically, fitting the mass difference Ec — Ac) play a 
minor role. In the following numerical calculations we will neglect the exchange 
of heavy mesons. 

Finally one has to fix qyajck masses. The light quarks ones are fixed to 
0.34 GeV according to Ref.sESH. The heavy quark masses mg — rric and mb 
are adjusted to reproduce the experimental average mass M = {M + 3M*)/4 
between the 0~ and 1~ M = D ot B mesons by a variational calculation, where 
a trial wave function of type (f> oc exp(— ar^/2) is used (with a as a variational 
parameter). It has been checked that the error never exceeds a few MeV 
with respect to the exact value. The variational approximation is retained for 
consistency with the treatment of 3-, 4- and 6-body systems discussed below. 
This leads to rric = 1.35 GeV and mb = 4.66 GeV for the potential Ci and 
mc = 1.87 GeV and rrif, = 5.259 GeV for Cj. 

Before proceeding further, let us briefly discuss the calculation of baryons 
masses in the model of Glozman et al. The explicit form of the Hamiltonian 
integrated in the spin-flavour space is : 

- 2 (go/g)' K,' for N 
H = Ho + -r^< ' , (8) 




with 



2 (go/5) K,' for A 



2m 2 



where g^/47r — 0.67 (which leads to the usual strength g7rJVAf/47r ~ 14 for the 
Yukawa tail of the nucleon-nucleon potential) and [ga/gY — 1-8- We have 
performed variational estimates with a wave function (j) cx exp(— a(p^ +A^)/2), 
where p = ?2 ~ ?3i ^ — (2?i — 1^2 ^ f3)A/3: our results agree with the more 
elaborated Faddeev calculations of Ref. E3. 

When the meson-exchange terms are switched off, the N and A ground 
states are degenerate at 1.63 GeV. Introducing the coupling the nucleon mass 
drops stronger leading to a reasonable splitting (w 0.3 GeV). It is interesting 
to notice that in this model one has an attraction both for the 5 = 0, / = 



and also, albeit smaller, for the 5 = 1, / = 1 diquarks. Then both the 
nucleon and A masses decrease when the interaction is turned on. The diquarks 
S = 0, I = 1 and S' = 1, / = on the contrary are repulsive configurations. 
The situation is thus different from the chromomagnetic case where one has a 
smaller attraction for the S = 0, I — diquark and repulsion for the S = 1, 
1=1 one. 

We have also calculated the ground state of cqq baryons using a trial wave 
function (j) cx exp(-(ap2 + i3X'^)/2) and found = 2.32 GeV and = S* = 
2.48JjeV, close to the experimental values and consistent with the findings of 
Rcf.LJ, although the Hamiltonian, its treatment, and the input parameters are 
somewhat different there. 



3 Evaluation of multiquarks masses 

Due to arguments at the beginning of this contribution, here we discuss tetra- 
quarks containing heavy flavours, i.e. QQqq, studying the most favourable 
configuration: 33, S' = 1, / = 0. This means that QQ is in a 3 colour state 
and qq in a S colour state. The mixing with 66 is neglected because one expects 
this to play a negligible role in deeply-bound heavy systemsB. Then the Pauli 
principle requires S12 = 1 for QQ, and 5*34 = 0, I34 = for qq (which we 
have seen to be the diquark with the largest binding), if the relative angular 
momenta are zero for both subsystems. This gives a state of total spin 5=1 
and isospin / = 0. 

The tetraquark Hamiltonian integrated in the colour-spin-flavour space, 
and incorporating the approximations discussed in the former section, reduces 
to 

H = 2(m + ruQ) + ^ + ^ + ^^±^p^ + ^ (10) 
niQ m imniQ ^ — ' 



where 



V12 = \ ( + hri2 + c 
2 V ^12 

= i ^--^ + 6ry , i = 1 or2, j ==3or4, (11) 

F34 = ^ ( -— + &r34 + C ) + - Vr, - 2Vr,'. 

2 V ^34 / 

The momenta p^,, etc., are conjugate to the relative distances x = ri — r2, 
y fs — r4, and z = (fi + f2 ~ ?3 ~ ?4) /v^- The wave function is parameterized 
as 

oc exp[-(aa;^ +/3y^ +7X^)/2], (12) 



and the minimization with respect to a, [3 and 7 shows that both the ccqq and 
hhqq systems are bound whatever is the potential, (Ci) or (C2), provided meson 
exchange is incorporated. This is in contradistinction to previous studies based 
on conventional models where the flavour-independent confining potential is 
supplemented by one gluon exchange. For example the authors of Rem found 
that the ccqq state is about 20MeV above threshold (while_the hbqq is bound 
of 135 MeV). A similar situation is also obtained in Ref. E3, where studying 
the four quarks states in a diquark model only the bbqq is found to be bound 
(by 50MeV, while ccqq is 30 MeV above the threshold). 

Quantitatively in our model we find that using the confining potential Ci 
plus meson exchange the double charmed tetraquark is bound by 185 MeV 
and the double bottom one by 226 MeV. Using the potential C2 the results 
are nearly two times larger (332 and 497 MeV respectively) and also much 
more different from each other. The reason is that (C2) contains a Coulomb 
part which binds more, heavier is the system, leading thus to a larger separa- 
tion among levels as well. This is related to the fact that the potential (C2) 
has been fitted to reproduce the J/^E* and the T meson masses (anyway it 
also gives ovecall good results both for other mesons and baryons) while, by 
construction EJ'Ej, the potential (Ci) was designed and fitted to light baryons 
only. Altogether, our results shows that the prediction of binding for the ccqq 
is substantially independent on the choice of the confining potential. 

Considering this result one could rise the question if in the model of Gloz- 
man et al. a proliferation of multiquark systems appears. We have therefore 
tried to investigate QQqqqq and q^ systems as well. 

As a general procedure, for a given multiquark system, one searches for 
the spin-isospin wave functions corresponding to a colour singlet, then selects 
the most favourable configuration. The contribution of global spin-fiavour- 
averaged interaction is reduced to the calculation of the matrix elements of the 
two body operator {ai-cTj)[fi-fj); this is accomplished by using Clebsh-Gordan 
coefficients of the permutation group according to Ref.sEao. 

Let us begin with QQqqqq. In this case the most favourable configuration 
is the one where the light quark subsystem has = 1, / = 0, which leads to 
the spin-flavour-averaged interaction 

{V) = IW. - 2/31/,, - 4/3(5o/g)V„'. (13) 

In the numerical calculation we have used the variational Gaussian wave 
function 

cx exp[-(Q;a::^ -|- Py"^ + j{u'^ + + u)^)] (14) 

where appear the Jacobi variables x = fi — r2, y = Rq — Rq, in terms of 
Rq = (?! -I- r2)/2 and Rq = (fa -I- f*4 -I- Fs -I- f*6)/4, u = (rs -f r4 - rs - r6)/\/2. 



V = (?3 ^ ?4 + ?5 ~ ?6)/'v/2, w = (13 — r4 — Fs + 1(,)/\/2^ where indices 1, 2 refer 
to the heavy quarks and 3, 4, 5, 6 to the hght ones. 

The result of the numerical calculation is that this potential is largely 
insufficient to bind both the ccqqqq (e.g. of 500 MeV with the potential Ci) 
and the bbqqqq systems. Incidentally, at this workshop results about heavy 
hexaquarks were presented also by Lichtenberg et al. 123, in their diquark model. 
Albeit they do not treat directly the case ccqqqq, they find that the system 
csqqqq is unbound, while the one containing a beauty quark instead of the 
charm is bound. 

For the sake of completeness we have also made a study of the q^ system. 
The most favourable configuration is 5 = 1, / = 0, leading to 

{V) = nV^ - 5/3K, - 10/3(go/g) (15) 

Our result shows that, also in this case, the system is largely insufficiently 
bound (by almost 1.5GeV with the potential Ci) for being under the two 
baryons threshold. 

4 Conclusions 

In conclusion, considering the success of the Glozman et al. model in describing 
the baryon spectrum, we have searched for new possible tests of this model. 
Unluckily the model does not permit an analysis of mesons (light-quark mesons 
at least), which are interpreted as "quasiparticles" related to the breaking of 
flavour SU{^). Further tests of it would then involve multiquarks systems. 
We have thus studied four and six-quark systems involving two heavy quarks, 
showing that the Glozman et al. model leads to predictions for the ccqq which 
differ from the ones of more conventional models. In fact this system is found to 
be strongly bound. The search of such a resonance in forthcoming experiments 
will therefore be a good test for understanding the dynamics which produces 
the hadron spectrum. 

We have also found that the bbqq state is bound, but this result is obtained 
also in other models and the possibility to observe this resonance is relegated 
to a more remote future. 

Finally we have found that the six quark system with two heavy quarks is 
unbound, also this result does not differ from the one of other more conven- 
tional models. However, together with the outcome that also the q^ system is 
unbound, it is a useful indication that in the framework of the model under 
investigation there is no proliferation of multiquark bound states. 



Acknowledgements 



We thank D.B. Lichtenberg for cnlighting discussions and the organizers for 
the stimulating atmosphere of this workshop. 



1. R.L. Jaffe, Phys. Rev. D 15, 267 (1977); Phys. Rev. D 17, 1444 (1978). 

2. R.L. Jaffe, Phys. Rev. Lett. 38, 195 (1977). 

3. A. T. M. Aerts et al, Phys. Rev. D 17, 260 (1977). 

4. J. Weinstein and N. Isgur, Phys. Rev. D 27, 588 (1983), Phys. Rev. D 
41, 2236 (1990). 

5. S.Zouzou, B. Silvestre-Brac, C. Gignoux and J.-M. Richard, Z. Phys. C 
30, 457 (1986). 

6. B. Silvestre Brae and C. Semay, Z. Phys. C 59, 457 (1993)); Z. Phys. 
C 61, 271 (1994). 

7. D.M. Brink and Fl. Stancu, Phys. Rev. D 49, 4665 (1994). 

8. J. Carlson and V.R. Pandharipande, Phys. Rev. D 43, 1652 (1991). 

9. J.P. Ader, J.-M. Richard and P. Taxil, Phys. Rev. D 25, 2370 (1982). 

10. G. Karl, Int. J. of Mod. Phys. El, 491 (1992); Nucl. Phys. A 558, 
113c (1993); N.A. Tornqvist, Proc. of "Int. Europh. Conf. on High 
Energy Phys.", Brussels (Belgium), edit. J. Lemonne et al., 84 (1995); 
G. Landsberg, P/iys. of Atom. Nucl. 57,42 ( 1994). 

11. M. Genovese, Nuovo Cimento A 107, 1249 (1994) and references therein. 

12. Particle Data Group, Phys. Rev. D 54, 1 (1996). 

13. M.A. Moinester, Z. Phys. A 355, 349 (1996). 

14. D.M. Kaplan, Proc. Int. Workshop "Production and Decay of Hyperons, 
Charm and Beauty Hadrons", Strasbourg (France) September 5-8, 1995. 

15. COMPASS Collaboration (G. Baum et al.), CERN-SPSLC-96-14, 
March 1996. 

16. S. Pepin, Fl. Stancu, M. Genovese and J.-M.Richard, |hep-ph 9609348| , 
ISN-96.99, to be published in Phys. Lett. B. 

17. S. Pepin, Fl. Stancu, M. Genovese, J.-M.Richard and S. Zouzou, work 
in progress. 

18. L.Ya. Glozman and D.O. Riska, Phys. Rep. 268, 263 (1996). 

19. L.Ya. Glozman, Z. Papp and W. Plessas, Phys. Lett. B 381, 311 (1996). 

20. L. A. Copley, N. Isgur and G. Karl, Phys. Rev. D 20, 768 (1979); S. 
Capstick and N. Isgur, Phys. Rev. D 34, 2809 (1986); C. S. Kalman 
and D. Pfeffer, Phys. Rev. D 27, 1648 (1983). 

21. N. Tornqvist, Phys. Rev. Lett. 67, 556 (1991); Z. Phys. C 61, 525 
(1994). 



22. A.V. Manohar and M.B. Wise, Nucl. Phys. B 399, 17 (1993). 

23. T.E.O. Ericson and G. Karl, Phys. Lett. B 309, 426 (1993). 

24. M. Weyrauch and H.J. Weber, Phys. Lett. B 171, 13 (1986); H.J. Weber 
and H.T. Williams, Phys. Lett. B 205, 118 (1988). 

25. L.Ya. Glozman and D.O. Riska, Nucl. Phys. A 603, 326 (1996). 

26. A. Buchman et aL, Nucl. Phys. A 569, 661 (1994). 

27. D. B. Lichtenberg, R. Roncaglia and E. Predazzi, IUHET-344, published 
in these proceedings. 

28. Fl. Stancu, Group Theory in Subnuclear Physics, Oxford University 
Press, 1996, chapter 4; 

29. S. Pepin and Fl. Stancu Preprint ULG-PNT-96-1-J. 

30. D. B. Lichtenberg, R. Roncagha and E. Predazzi, DFTT 65/96, published 
in these proceedings. 



